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^ ; Abstract 

' In models involving new TeV-scale Z' gauge bosons, the new U{1)' symme- 

. try often prevents the generation of Majorana masses needed for a conventional 

neutrino seesaw, leading to three superweakly interacting "right-handed" neu- 
trinos lyR, the Dirac partners of the ordinary neutrinos. These can be produced 
prior to big bang nucleosynthesis by the Z' interactions, leading to a faster 
expansion rate and too much ^He. We quantify the constraints on the Z' prop- 
erties from nucleosynthesis for Z' couplings motivated by a class of Eq models 
. parametrized by an angle 9 eg- The rate for the annihilation of three approxi- 

mately massless right-handed neutrinos into other particle pairs through the Z' 
channel is calculated. The decoupling temperature, which is higher than that 
of ordinary left-handed neutrinos due to the large Z' mass, is evaluated, and 
5^ I the equivalent number of new doublet neutrinos ANu is obtained numerically 

as a function of the Z' mass and couplings for a variety of assumptions concern- 
ing the Z — Z' mixing angle and the quark-hadron transition temperature Tc- 
Except near the values of 9 eg for which the Z' decouples from the right-handed 
neutrinos, the Z' mass and mixing constraints from nucleosynthesis are much 
more stringent than the existing laboratory limits from searches for direct pro- 
duction or from precision electroweak data, and are comparable to the ranges 
that may ultimately be probed at proposed colliders. For the case Tc = 150 
MeV with the theoretically favored range of Z — Z' mixings, AA'^^ ^ 0.3 for 
Mz' ^ 4.3 TeV for any value of 9e6- Larger mixing or larger Tc often lead to 
unacceptably large AA^^, except near the decoupling limit. 
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1 Introduction 



Additional heavy Z' gauge bosons P!J are predicted in many superstring and grand 
unified jH] theories, and also in models of dynamical symmetry breaking |3]. If present 
at a scale of a TeV or so they could provide a solution to the /i problem j^j and other 
problems of the minimal supersymmetric standard model (MSSM) f^. Current limits 
from collider [HIH] and precision 9J experiments are model dependent, but generally 
imply that Mz> > (500 — 800) GeV and that the Z — Z' mixing angle is smaller than 
a few xlO^"^. There are even hints of deviations in atomic parity violation JHl^- and 
the NuTeV experiment ^21, which could be an early indication of a Z' JH]- A Z' 
lighter than a TeV or so should be observable at Run II at the Tevatron. Future 
colliders should be able to observe a Z' with mass up to around 5 TeV and perform 
diagnostics on the couplings up to a few TeV [14 . 

An electroweak or TeV-scale Z' would have important implications for theories 
of neutrino mass. If the right-handed neutrinos carry a non-zero 11(1)' charge, then the 
U{1)' symmetry forbids them from obtaining a Majorana mass much larger than the 
'-breaking scale, and in particular would forbid a conventional neutrino seesaw 
model In this case, it might still be possible to generate small Majorana masses 
for the ordinary (active) neutrinos by some sort of TeV-scale seesaw mechanism in 
which there are additional mass suppressions ^HI- However, another possibility is 
that there are no Majorana mass terms, and that the neutrinos have Dirac masses 
which are small for some reason, such as higher dimensional operators ^7] or volume 
suppressions in theories with large extra dimensions ^Hl- In this case, the model 
would contain three additional right-handed partners of the ordinary neutrinos, which 
would be almost massless. Such light Dirac neutrinos (i.e., with mass less than an 
eV or so) in the standard model or MSSM are essentially sterile, except for the tiny 
effects associated with their masses and Higgs couplings, which are much too small 
to produce them in significant numbers prior to nucleosynthesis or in a supernova. 
However, the superweak interactions of these states due to their coupling to a heavy Z' 
(or a heavy W in the SU{2)l x SU{2)r x U{1) extension of the standard model [TU] ) 
might be sufficient to create them in large numbers in the early universe pUl 1^ 122) 
or in a supernova j2S]- In this paper, we consider the constraints following from big 
bang nucleosynthesis on Z' properties in a class of E'e-motivated models. 

It is well known that any new relativistic particle species that were present 
when the temperature T was a few MeV would increase the expansion rate, leading 
to an earlier freeze-out of the neutron to proton ratio and therefore to a higher 
'^He abundance [211 12H] • Their contribution is usually parametrized by the number 
AA^,^ of additional neutrinos with full-strength weak interactions that would yield 
the same contribution to the energy density. The primordial '^He abundance is still 
rather uncertain, but typical estimates of the upper limit on AA^^, are in the range^ 

^The interpretation of these results is controversial. For recent discussion, see 
^The limit can be weakened by invoking an excess of Vf, with respect to Pg, which lowers the n/p 
ratio. 
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ANy < (0.3 — 1) fIE\ 1^ . Of course, the Z-width does not allow more than 3 
light active neutrinos j^ZI; so AA^",^ should be interpreted as an effective parameter 
describing degrees of freedom that do not couple with full strength to the Z. 

In 1979, Steigman, Olive, and Schramm [2()\ 121] described the implications of 
a superweakly interacting light particle, such as a right-handed neutrino coupling to 
a heavy Z'. Because of their superweak interactions, such particles decoupled earlier 
than ordinary neutrinos. As the temperature dropped further, massive particles such 
as quarks, pions, and muons subsequently annihilated, reheating the ordinary neu- 
trinos and other particles in equilibrium, but not the superweak particles. One must 
also take into account the transition between the quark-gluon phase and the hadron 
phase. 

A simple estimate of the decoupling temperature is obtained as follows [201 In- 
ordinary neutrinos have cross-sections aw oc G^T"^, where Gw is the Fermi constant, 
and interaction rates 

TwiT) = n (awv) cc GI,T\ (1) 

where n is the density of target particles. The Hubble expansion parameter varies 
as if oc T"^ /Mp, where Mp is the Planck scale, so the decoupling temperature at 
which r is equal to H becomes 

T, oc {Gl,Mp)-''\ (2) 

Putting in the coefficients, Tct{vi) k, 1 MeV for the ordinary neutrinos^. Similarly, a 
superweakly interacting particle such as a right-handed neutrino with a cross-section 
asw GliyT^, would decouple at 

T.M - l^) Uvl). (3) 

If in the specific model, the effective superweak coupling constant Gsw is propor- 
tional to iW^sw, where Mgy/ is the mass of a superweak gauge boson, the decoupling 
temperature can be written as 

where My/ is the W mass. It is then straightforward to calculate the dilution by the 
subsequent quark-hadron transition and the annihilations of heavy particles, and the 
corresponding AA^;^ from the superweak particles. 

Of course, the estimate in (jU is very rough. In particular, the detailed cou- 
plings of the Z' to the and to all of the other relevant particles must be considered 

■^More detailed studies obtain Td{i^L) ^ 3 MeV. We will obtain Td{vR) by an explicit calcu- 
lation, so the difference is irrelevant for our purposes. 
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for a precise estimate^ In this paper, we do this for a class of Z' models with couplings 
motivated by grand unification (The full structure of is not required.) 

We define the UiX)' model in Section |21 The implications of superweakly coupled 
particles for nucleosynthesis and the uncertainties from the quark-hadron transition 
temperature are summarized in Section El Section E] deals with the calculation of 
the decoupling temperature. We present our results and numerical analysis for and 
lS.Ny for three right-handed neutrinos as a function of the Z' mass and couplings for 
various assumptions concerning theZ — Z' mixing and in Sectional The discussion 
and conclusion follows in Section IHl 

2 Z' in E'e-motivated models 

A general model with an extra Z' is characterized by the Z' mass; the Z — Z' mixing 
angle; the gauge coupling; the chiral charges for all of the fermions and 

scalars, which in general may be family non-universal, leading to flavor changing 
neutral currents [33]; and an additional parameter associated with mixing between 
the Z and Z' kinetic terms Furthermore, most concrete Z' models involve 

additional particles with exotic standard model quantum numbers, which are required 
to prevent anomalies. It is difficult to work with the most general case, so many 
studies make use of the charges and exotic particle content associated with the 

E^ model, as an example of a consistent anomaly-free construction^. Explicit string 
constructions jHS] often lead to other patterns of couplings and exotics, but these are 
very model dependent. 

E^ actually yields two additional f/(l)' factors when broken to the standard 
model (or to SU{5)), i.e., 

Ee SO(IO) X U{1)^ SU{5) x f/(l)^ x f/(l)^. (5) 

It is usually assumed that only one linear combination survives to low energies, 
parametrized by a mixing angle Beg- The resultant U{1)' charge is then^ 

Q = Q^cos6e6 + Q^sm6E6. (6) 

A special case that is often considered is U{l)r,, which corresponds to 9e6 = 2n — 
ta.n'^ = l.Tlvr. We list the charges of U{1)^ and U{1)^ that we need in Table 
n The quantum numbers of the associated exotic particles are given in j^. It is 

""^Detailed calculations were carried out in "SF for the 77 model (see Section EJ, in [221 for more 
general Eq models, and in 30 for a model with generators T^r and B — L. However, these studies 
considered only vrPe. (e+e^, vlPl)- In the present paper we include the interactions with all of 
the particles in equilibrium at a given temperature. This leads to a lower Td{i'R) and more stringent 
limits. Constraints on extended technicolor models were considered in 'Jll . 

^The full structure of Eq grand unification is not required, and in fact the Eq Yukawa coupling 
relations must not be respected in order to prevent rapid proton decay j32| . 

^ We ignore the possibility of kinetic mixing [31] . 
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Table 1: The (family-universal) charges of the U{1)^ and the U{1)^. 



Fields 


Qx 




Ul 


-l/2v^ 


1/2^6 


Ur 


l/2v^ 


-l/2v^ 


dL 


-l/2v^ 


1/2^6 


dR 


-3/2v^ 


-l/2x/6 




3/2v^ 


l/2\/6 


dR 


l/2v^ 


-l/2v^ 




3/2710 


l/2v^ 




5/2v^ 


-l/2v^ 



conventional to choose 9e& to be in the range (0, tt), since the charges merely change 
sign for Oe^ — > Oe& + tt. With this convention one must allow both positive and 
negative values for the Z — Z' mixing angle 5. In this paper, we find it convenient 
to choose a different convention in which 9e& varies from to 27r, but for which 
5 < 0. That is, the range — tt corresponds to the models with negative mixing, 
while vr — 27r corresponds to positive mixing. The vr charge is nonzero, precluding 
an ordinary seesaw, except for 9eg ^ 0.427r and 1.427r. We will always assume that 
the neutrinos are Dirac and that the three right-handed neutrinos are therefore very 
light. (In fact, the non-zero Dirac masses play no role in the analysis.) There could 
be additional sterile states, such as the S'O(10)-singlet states occurring in the 27-plet 
of Eq. If these involve nearly- massless fermions they could also contribute to the 
expansion rate prior to nucleosynthesis. We assume that these additional neutralinos 
acquire electroweak scale masses from the gauge symmetry breaking pT^ . 

Let Z and Z' represent the Standard Model and U{1)' gauge bosons, respec- 
tively, and Zi 2 the mass eigenstate bosons, related by 

^ Zi\ ( cos 5 — sin (5 \ / Z 
^ ^2 J y sin5 cos 5 J \ ^' 

where 6 is the Z — Z' mixing angle. As stated in the Introduction, the limits on 
~ Mz' depend on 9 eg and also on the masses of any exotics and superpartners to 
which the Z' couples, but are typically in the range Mz' > (500—800) GeV. The limits 
on 6 are correlated with those for Mz' and are asymmetric under 6 —6. However, 
for Mz' ~ 1 TeV the constraints are less sensitive to 9 eg and are approximately 
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symmetric, with \6\ < 0.002 giving a reasonable approximation for all 6e6- For larger 
Mz' there are two theoretical constraints on the mixing, corresponding to equations 
(6) and (5) of fSUj. The first is a theoretical relation between the mass and mixing. 



where Qz = y 9i + 92 and g'^ is the U{1)' gauge coupling constant. The value of g'^ 
depends on the embedding and breaking of the underlying theory. We will choose 
g'z = ^J^gz sm9w, which corresponds to a unification of g'^ with the other gauge 
couplings for the exotic particle quantum numbers of supersymmetric Eg. In (jHl) C 
depends on the charges of the scalar fields which lead to the mixing (see Table III 
of [SS])- However, for the typical cases in which the mixing is induced by scalars in 
an Eq 27 or 27-plet, it is a reasonable approximation to take — 1 < C < 1 for all 
6e6- (One can have a slightly more restrictive range for some Beg-) The assumption 
|C| < 1 corresponds to\6\ < O.OOSl/Mf^, where M^^ is in TeV. The second theoretical 
constraint is the requirement that the mixing should not change the mass of the lighter 
Z more than is allowed by the data. It is equivalent to 

M: 



\S\-^Po-l (9) 

where Mz^ = Mz, and the po parameter, defined precisely in [Slj, should be ex- 
actly 1 in the standard model. The precision data imply po < 1.001. Hence, 
\S\ < 0.0029/M^2, where Mz2 is again in TeV. We will consider the following cases: 

(AO) 6=0 (no mixing) 

(Al) \6\ < 0.0051/M|^ (theoretical mass — mixing relation) 

(A2) \6\ < 0.0029/M^2 (po constraint) 

(A3) \6\ = 0.002 (maximal mixing allowed for Mz^ ~ 1 TeV). (10) 

Al is more stringent than A2 and A3 in the large mass range, so we will mainly focus 
on AO and Al. 

The lagrangian for the massive neutral current coupling to fermion / is j^ni 

-^int = 9zQzifL)fLYfLZ^,+ gzQz{fR)fRYfRZi, 

+9'zQifL)fLYfLZ'^ + 9'zQifR)fRrfRZ', (11) 

where 

Qz{fL) = Tf-qfSm^9w, 

Qz{fn) = -qfSmHw, (12) 

and Q{fL,R) is given by ©. The annihilation cross-section through Z' has both (light) 
Zi and (heavy) Z2 contributions unless 6 = and is calculated in Section |3J 
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3 Nucleosynthesis 

As described in the Introduction, the observed ^He abundance constrains the energy 
density at the time of Big Bang Nucleosynthesis |21I, with most recent estimates [2^1 
EE] of the number of equivalent new active neutrino types in the range AN,y < (0.3—1). 
The contribution of new relativistic species can be written 



where qb and gp are degrees of freedom of new bosons (B) and new fermions (F), 
respectively, Tb^f are their effective temperatures, and Tbbn ~ 1 MeV is the tem- 
perature at the time of the freeze-out of the neutron to proton ratio. In particular, 
the contribution of three types of right-handed neutrinos is 



where Tdiyn) is the decoupling temperature of the right-handed neutrinos. g{T) 
is the effective number of degrees of freedom at temperature T. Neglecting finite 
mass corrections, it is given by 5's(T) -|- ^gF(T), where gB,F(T) are the number of 
bosonic and fermionic relativistic degrees of freedom in equilibrium at temperature 
T [201 EI]- In particular, g(TBBN) = 43/4 from the three active neutrinos, e^, and 
7, and g{T) increases (in this approximation) as a series of step functions at higher 
temperature as more particles are in equilibrium. The second equality in (jl4|) comes 
from entropy conservation pUl in the heavy particle decouplings and quark-hadron 
transition subsequent to the uji decoupling. Therefore, the are not included in 
our definition of g{T). (They will be included in the expansion rate formula prior to 
decoupling.) 

In calculating g{T) one must also take into account the QCD phase transition 
at temperature Tc. Above Tc the u and d (and possibly s) quarks and the gluons 
were the relevant hadronic degrees of freedom, while below Tc they are replaced by 
pions [201 1^ ■ The value of is poorly known, but is usually estimated to be in the 
range (150 — 400) MeV j2Hl- This range is estimated in quark and hadron potential 
models as the temperature above which hadrons start to overlap (lower end) or as 
the temperature below which the quark gas in no longer ideal (upper end). A related 
uncertainty is whether to use current or constituent quark masses. At very high 
temperatures the quarks can be considered as asymptotically free and current masses 
are appropriate, while around Tc constituent effects become important''. The range 
of estimates for Tc is essentially unchanged if one simply fixes the quark masses at 
either value [5^ . 

^One can alternatively argue that the current masses are appropriate above a temperature Tchirai, 
above which chiral symmetry is restored, and constitutent masses below Tchirai- One would expect 
Tc and Tchirai to be comparable, but their precise relation is uncertain. 




(13) 




(14) 
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Figure H shows the exphcit values of g{T) from the more detailed analysis of 
Ref. which includes finite mass and other corrections, and uses the two values 
Tc = 150 MeV and 400 MeV. We will also use these values for our numerical analysis. 
The sharp increase in g{T) above Tc (because of the large number of quark and gluon 
degrees of freedom) is extremely important for relaxing the constraints on the Z' 
mass. 

The QCD phase transition does not occur instantaneously or at one tem- 
perature but rather smoothly (meaning both quarks and hadrons exist at the same 
temperature) for a period of time around Tc, as illustrated by the smooth curves in 
Figure ^ Risking a small inconsistency, we approximate our calculation of the inter- 
action rate by simply switching from quarks to hadrons for temperatures below Tc. 
We will take the values Tc = 150 and 400 MeV to illustrate the range of hadronic un- 
certainties. Above Tc, the interaction rate depends in principle on the quark masses, 
especially for low Tc. However, we have found in practice that the results are almost 
identical for constituent and current masses, so we will mainly display them for the 
constituent case (both will be shown for the rj model). 

The calculation of the right-handed neutrino decoupling temperature, Td^i'R) 
in terms of the Z' parameters is discussed in the next section. 



4 The expansion and interaction rates 

A particle is decoupled from the background when its interaction rate drops below 
the expansion rate of the universe. In this section, we present the the cosmological 
expansion rate H{T) along with the explicit form of the interaction rate r(T) for 
TTjii^R annihilating into all open channels^, and estimate the decoupling temperature 
Td of a right-handed neutrino by T{Td) ~ H(Td). 

The Hubble expansion parameter is given by 

^^^^^^S^G^^^I^Ggm^. (15) 

where Gn = Mp'^ is the Newton constant and p(T) is the energy density. We define 
g'{T) = g{T) + ^, where the 21/4 reflects the 3 massless right-handed neutrinos. 

The cross-section (ri{s) = crijTRVji fifi) for a massless right-handed neutrino 
pair to annihilate into a fermion pair through the Z'-channel is 



®As long as equilibrium is maintained, the vn annihilation and production rates are the same. It 
is more convenient to estimate the annihilation rate oiTT^vn into massive particles, because the final 
state mass effects are easily incorporated in the cross section formulae, whereas for the production 
rate one must explicitly consider the suppressed number density for the massive particles. 
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where (for s < M|^, M|J 



- gz9zQ[vR)QzUix)\^T-2 T72 I' (17) 



where X = L or i?, Pi = ^1 — Amjjs is the relativistic velocity for the final 

particles, and A^^ is the color factor of particle fi. 

In the limit of no-mixing {6 = 0) and massless final particles (/3j = 1), the 
cross-section simplifies to 



-J2 



2 



^^^'^ ^ [j^] QM' [Qif'^LY + Qif, 



/2 

consistent with the earlier estimate asw G'g^rT'^ with Gsw and T oc ^/s. 

For temperatures less than the quark-hadron transition temperature = 
150 — 400 MeV, we replace the quark degrees of freedom with hadrons. The only 
relevant annihilation channels are into charged pions. We approximate the cross- 
section of TT^i^R annihilating into n~^7r~ by using the p dominance model 



a^{s) ^ a{T7Eun tt+tt-) = ^\FAs)\' (G^ + + G^j, + G^)' (19) 

which is basically obtained by using Q^fn) = Q{uL) + Q{dL) and Qzifn) = Qz{ul) + 
Qzidi) for G^j;^ and likewise for G^^. The pion form factor^ is 

2 

TfL 

= ^, (20) 



s — ml + impTp ' 

.e 

The interaction rate per uji is 



with mp = 771 MeV and Tp = 149 MeV. 



r(T) = Y: r.(T) =Y^!^ (av{T7Eun - vr+Tr")) , (21) 

where riy^ is the number density of a single flavor of massless right-handed neutrinos 
plus antineutrinos, = 2 is the number of degrees of freedom, and {av) is the 
thermal average of the cross-section times velocity. 

We use the same masses (Table |2I) used in the calculation [211 EHI of g{T) 
in Figure n except for the value mf, = 4200 MeV of the h quark current mass |27j . 
We include the contributions of all particles up to the h quarks. The contributions 

^More complicated form factors are known to fit the experimental data better but H2U|) is 
adequate for our purposes. 
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Table 2: The masses (in MeV) used for the numerical analysis. 



Quarks 


Current (Constituent) masses 


Others 


Masses 


u 


4.2 (340) 


V 





d 


7.5 (340) 


e 


0.511 


s 


150 (540) 




105 


c 


1150 (1500) 


r 


1800 


b 


4200 (4500) 


TT 


137 



from the top quark and heavy particles from new physics, such as squarks, sleptons, 
and exotics would only be relevant when the decoupling temperature is close to the 
electroweak scale or higher. This only occurs when is extremely close to the values 
for which the z/^ decouples from the Z' . 

For a massless right-handed neutrino pair colliding with 4-momenta = (p, p) 
and = (fc, k) with relative angle 6', the interaction rate per neutrino is 

_ 9un r^2,^r,2„f'^^^n (l-COS( 



fp'dpf edkl^dcose ^ ^^^^ ^ ^^ a,{s), (22) 



where fuik) = (e*^/^ + 1) ^ is the Fermi-Dirac distribution with 

d^k 3 



/a k -i 
(^/^(^) = 2 ■ (23) 



vm = P ■ k/pk = 1 — cos^ is the M0ller velocity, and s = 2pk{l — cos 6') is the square 
of the center-of-mass energy. 

A root-finding method was used to calculate the decoupling temperature, for 
which H = r. A several percent error was allowed in the numerical result to calculate 
the roots efficiently. Finite temperature effects, such as changes in the phase space due 
to interactions with the thermal bath, can increase the ordinary neutrino decoupling 
temperature by several percent [121 • Analogous effects for the are too small to 
significantly affect our results. 



5 Numerical results 

In this section, we present the numerical results from the calculation. The marked 
points in Figures |2E are the results of the actual calculation, while the curves inter- 
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polate. 

Figures El and El show how the right-handed neutrino decouphng temperature 
Trf and the equivalent number of extra neutrino species AA*",^ change with for 
9e& = 27r — tan~^ ~ l.Tlvr (the rj model) for constituent and current masses, 
respectively, for = 150 and 400 MeV and the various assumptions concerning the 
Z — Z' mixing listed in The no- mixing curves (AO) exhibit an approximate 

Td ~ {MzjMzf^ dependence, in agreement with the simple estimate in the Intro- 
duction [2ni 121] • This is to be roughly expected because of the M^^ dependence of 
the cross section for no mixing, but is not exact because additional channels which 
affect both the expansion and interaction rates open up at higher temperatures. The 
no- mixing curves in Figures El and El are reasonably described by Q for Tfi{i'L) ~ 
3 MeV for the t] model, but the coefficients in front of (M^j/^z)^^^ strongly 
model dependent, as is apparent in Figures OSl is usually lower in the cases 
involving Z — Z' mixing, because the Z annihilation channel yields a contribution 
proportional to 5^ even for infinite M^j- That is why the (theoretically unrealistic) 
curves A3 for fixed \5\ = 0.002 are asymptotically fiat for large M^j. Case Al, in 
which \6\ ~ 0.0051/M4, also has Ta ~ {MzjMzY^^, though with a smaller coeffi- 
cient than for no mixing^°, while A2, with \5\ = 0.0029/M^2, has Td ~ (MzjMz)^^^. 
For case Al, Td is asymmetric under 6 —6 for all Mz2, as is apparent from (fTUI) 
and p7p. The difference vanishes asymptotically for A2 and A3, but even for Mz2 = 
5 TeV there is still a difference, especially for A2. 

The decoupling temperature is slightly lower for = 400 MeV than for 150 
MeV, provided it is in the range for which the two curves in Figure Q differ. Both 
the expansion and annihilation rates are smaller for = 400 MeV, but the effect 
on the expansion rate is more important because of the gluonic degrees of freedom. 
Similarly, Td is smaller for current quark masses than for constituent masses, provided 
Td > Tc, because of the larger annihilation rate^^. 

The ANi, curves change rapidly when Td reaches the quark-hadron phase tran- 
sition temperature Tc, where g{T) changes significantly. That is why ANi, is so much 
larger for Tc = 400 MeV than for 150 MeV. For the no-mixing case, the difference is 
significant for ^ 4 TeV, and it persists to even higher masses for the mixing cases 
(and to infinite mass for maximal mixing). The only significant difference between 
the constituent and current quark masses is in the maximal mixing case with Tc = 
150 MeV. That is because Td is very close to Tc, and even a small change in Td leads 
to a significant change in g{T), as can be seen in Figure ^ 

It is apparent from Figures El and 01 that the rj model leads to a significant 
AA^,^ for all of the cases and parameter ranges considered. Even the very conservative 
constraint ANi^ < 1 implies Mz2 > 1-5 — 2.2 TeV for Tc = 150 MeV, or, limiting 
ourselves to the most realistic cases AO and Al, Mz2 > 1-5 — 1.9 TeV. For Tc = 400 

^°The coefficient is smaller for most but not all values of 9e6- 

"'^"'^The difference between current and constituent masses would be reduced if their effects in the 
annihilation rate were properly correlated with those in the expansion rate. However, as described in 
Sectional the effect on g{T) is small compared with the uncertainty from Tc, and will be neglected. 
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MeV one finds Mz^ > 3.3 -4 TeV for AO and Al, Mz^ > 5 TeV for A2 and no allowed 
values for A3. All of these are much more stringent than the direct laboratory limit 
of 620 GeV [7j or the indirect limits from precision electroweak data The more 
stringent limit AA'^^ < 0.3 is satisfied for cases AO and Al for Mz2 > 2.5 — 3.2 TeV 
for = 150 MeV, and Mz^ > 4.0 - 4.9 TeV for = 400 MeV. It is not satisfied 
for case A2 with Tc = 400 MeV until extremely high masses, and never for (fixed) 
maximal mixing unless one takes a mixing much smaller than the present accelerator 
limit {\6\ < 0.0024) 8J. 

Figures E] and IHl display the results for the class of Eq models parametrized by 
the angle Oeq defined in for constituent masses and = 150 MeV and 400 MeV, 
respectively. Each figure includes the no-mixing case and the mixing assumption Al 
defined in ()10p. which is the most stringent and realistic. The limits in the presence of 
Z — Z' mixing are asymmetric under S — > —S. This is represented in the right-handed 
graphs by taking 6 < but allowing 6e6 to run from to 27r, so that the (vr — 27r) 
range for 5 < is equivalent to (0 — tt) with 6 > 0. The top graphs display as 
a function of Oeq for fixed values M^, = 500, 1000, 1500, 2000, 2500, 3500, 4000, 
and 5000 GeV, with larger Mz2 corresponding to higher T^. The middle graphs show 
AN,^ as a function of 9e6 for the same values of Mz2, with larger corresponding 
to smaller ANi,. The bottom figures show the lower bounds on for AA"^ < 0.3, 
0.5, 1.0 and 1.2, with larger AA^^ corresponding to smaller M^j. 

It is seen that becomes very large and the limits essentially disappear as 
9e6 approaches 9e6 ~ 0.427r or 1.427r, for which decouples completely {Q{i^r) = 0), 
but the details depend on the new physics at the electroweak and higher scales (we 
only explicitly included particles up to the b quark). 6eq = l.Tlvr corresponds to the 
f] model with 5 < 0, while 9eq = O.Tlvr corresonds to 5 > 0. It is seen from the figures 
that AA",^ is larger for values of Oeq closer to (the x model), but are weaker near 
Oeq = 7r/2 (the model). 

From the figures it is apparent that requiring AN^, < 1 excludes much of the 
interesting parameter space for Tc = 150 MeV, except for large Z2 masses or regions 
very close to the Ufi decoupling angles ~ 0.427r and lA2n. In particular, the AA^^ < 1 
constraint is satisfied for all values of 6e6 for ^ 2.2 TeV if there is no mixing, 
with a slightly more stringent constraint ^2.4 TeV for mixing assumption Al. 
The corresponding limits for AA^^ < 0.3 are 3.8 and 4.3 TeV. The constraints 
for Tc = 400 MeV are even more stringent, essentially requiring i/^ decoupling or 
very large Z2 masses. One has AA^j^ < 1(0.3) for all 6e6 for cases AO and Al for 
Mz2 ^ 5.1(6.1) TeV. 

6 Discussion and Conclusion 

Many theories beyond the standard model predict the existence of additional Z' gauge 
bosons at the TeV scale. The associated f/(l)' gauge symmetry often prevents the 
large Majorana masses needed for an ordinary neutrino seesaw model. One possibility 
is that the neutrino masses are Dirac and small. In that case, there is a possibility 



12 



of producing the sterile "right-handed" neutrino partners via Z' interactions prior 
to nucleosynthesis [201121], leading to a faster expansion and additional ^He. 

We have studied the right-handed neutrino decoupling temperature in a 
class of ii^6-niotivated U{1)' models as a function of the Z' mass and couplings (de- 
termined by an angle 9e%) for a variety of assumptions concerning the Z — Z' mixing 
angle 5, the quark-hadron transition temperature Tc, and the nature (constituent or 
current) of the quark masses. We have taken all relevant channels (quark, gluon, 
lepton, and hadron) into account, not only in the expansion rate H{T) and entropy, 
but also in the rate T{T) for a massless right-handed neutrino pair to annihilate into 
a fermion or pion pair via the ordinary or heavy Z bosons. We therefore obtain a 
larger annihilation rate, and thus a lower decoupling temperature and more stringent 
constraints, than earlier calculations, which only included annihilation into e^e~ and 

From the decoupling temperature and entropy conservation as quarks and 
gluons are confined or as various heavy particle types decouple and annihilate, one 
can obtain the number of right-handed neutrinos at nucleosynthesis, expressed in 
terms of the equivalent number AA^;^ of new ordinary neutrino species, for various 
sets of model parameters Mz^, 6eq, and Tc. Most recent studies of the primordial 
abundances obtain upper limits on AA^",^ in the range (0.3-1) As can be seen 

in Figures 131131 this implies rather stringent constraints on the Z' parameters for most 
values of Oeq- For = 150 MeV, the constraint AN,, < 0.3(1) is satisfied for all 6e6 
for Mz^ ^ 3.8(2.2) TeV for noZ-Z' mixing, and for Mz^ ^ 4.3(2.4) TeV allowing the 
range of mixing angles 6 obtained approximately when one assumes that the scalar 
fields responsible for the mixing are contained in the 27 or 27-plet of Eq (case Al in 
(I10|) ). For Tc = 400 MeV the constraints are much stronger, Mz2 ^ 6.1(5.1) TeV for 
AA*",^ < 0.3(1). The strong dependence on Tc is due to the large increase in the number 
of degrees of freedom for temperatures ^ Tc (Figure (d}), so that the number density 
of is strongly diluted for T^ ^ Tc. The constraints are strongest for 9e6 close to 
or TT, i.e., near the x model, which corresponds to S'O(IO) SU{5) x U{1))^, and are 
very weak near the ip model corresponding to Eq S'O(IO) x f/(l)^, Oeq = vr/2. They 
disappear entirely at the values Beg = 0.42-7? and 1.427r, for which the z/^j decouple 
from the Z'. The often considered r] model, 9e6 = 27r— tan^^ = l.Tlvr (or O.Tlvr for 
—Zfj) is somewhere in between, with the constraints shown in more detail in Figures 
El and El 

Except near the z//j decoupling angles, the Z' mass and mixing constraints 
from nucleosynthesis are much more stringent than the existing laboratory limits 
from searches for direct production or from precision electroweak data, and are com- 
parable to the ranges that may ultimately be probed at proposed colliders. They are 
qualitatively similar to the limits from energy emission from Supernova 1987A |23j . 
but somewhat more stringent for AN,, < 0.3, and have entirely different theoretical 
and systematic uncertainties. 

There are several ways to evade the nucleosynthesis constraints on an extra 
Z'. One possibility is to generate small Majorana neutrino masses for the ordinary 
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neutrinos by invoking an extended seesaw model ^j, in which the extra sterile neu- 
trinos are typically at the TeV scale. Another possibility is that the decouple from 
the Z', in which case the constraints disappear. This can in fact occur naturally in 
classes of models in which one combination of the x ci-nd ip charges is broken at a 
large scale associated with an F and D-flat direction leaving a light Z' which 
decouples from the ur^^. Yet another possibility is to weaken the observational con- 
straint on ANi, by allowing a large excess^^ of z/g with respect to z/g. This would, 
however, require a somewhat fine-tuned cancellation between the effects of the ur 
and the z/g — Ue asymmetry. 

Similar constraints on the W and Z' properties in SU{2)l x SU{2)r x U{1) 
models JH] are under investigation [IE] . 
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Figure 1: The effective number of degrees of freedom as a function of temperature 
for the quark-hadron transition temperature Tc = 150 MeV and 400 MeV, from |39j . 
g{T) does not include contributions from the three right-handed neutrinos, which are 
added separately in the expansion rate formula. 
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Figure 2: The decoupling temperature (top) and the equivalent number of ex- 
tra neutrinos AA^,^ (bottom) for the t] model as a function of the Z2 mass for 
constituent quark masses, for a quark-hadron transition temperature = 150 MeV 
(circles) and 400 MeV (crosses). The left two figures are for the cases AO and A3 
defined in (fTUI) . i.e., the solid, dashed and dotted lines represent zero- mixing (5 = 0), 
and positive and negative maximal-mixing {6 = ±0.002), respectively. The Tc = 150 
MeV case has higher and lower AN,y for the same than Tc = 400 MeV. The 
right figures are for the intermediate mixing assumptions Al and A2. The solid and 
dash-dot curves are for the mass-mixing relations 6 = ±0.0051/M|2, while the dashed 
and dotted curves are for the po constraints 6 = ±0.0029/M^2. 
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Figure 3: Same as Figure El except that current quark masses are used. The upper 
graphs share most features with the constituent mass case except that Td can be 
shghtly lower when Td > Tc. The only significant change in AN^y is for the Tc = 150 
MeV maximal mixing case (see text). 
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Figure 4: (top) and AA^"^ (middle) for Mz^ = 500, 1000, 1500, 2000, 2500, 3500, 
4000, and 5000 GeV, for Tc = 150 MeV and constituent masses. Larger cor- 
responds to higher Td and smaller AN^. The graphs on the left are for no mix- 
ing (case AO in ()10|l ). while the right-hand graphs are for the mass-mixing rela- 
tion \6\ < O.OO5I/MI2 (case Al). The bottom graphs are Mz^ corresponding to 
AN^ = 0.3, 0.5, 1.0 and 1.2, with larger AA^j, corresponding to smaller Mz2- 
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Figure 5: Same as Figure El except Tc = 400 MeV. is slightly smaller (for > 
150 MeV) for fixed and Beg, while AA^,, and the bound on Mz^ for fixed AA^^ 
are increased. 
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